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1. Introduction. 



The construction of explicit solutions to partial differential equations by symmetry 
reduction dates back to the original work of Sophus Lie, [30]. The reduction of partial 
differential equations to ordinary differential equations was generalized by Clarkson and 
Kruskal, [6], in their direct method, which was later shown, [35], to be included in the 
older Bluman and Cole nonclassical symmetry reduction approach, [3]. A survey of these 
methods as of 1992 can be found in [33]. Meanwhile, Galaktionov, [12], introduced the 
method of nonlinear separation that reduces a partial differential equation to a system 
of ordinary differential equations, developed in further depth in [12, 13, 14, 15]. Similar 
ideas appear in the work of King, [26], and the "antireduction" methods introduced by 
Fushchych and Zhdanov, [10, 9]. Svirshchevskii, [49, 50, 51], made the important obser- 
vation that one could, in the one- dimensional case, characterize in terms of higher order 
(or generalized) symmetries those nonlinear ordinary differential operators that admit a 
given invariant subspace leading to nonlinear separation. 

In quantum mechanics, linear differential operators with invariant subspaces form the 
foundation of the theory of quasi-exactly solvable (QES) quantum models as initiated by 
Turbiner, Shifman, Ushveridze, and collaborators, [41, 42, 43, 52]. The basic idea is that a 
Hamiltonian operator which leaves a finite-dimensional subspace of functions invariant can 
be restricted to this subspace, resulting in an eigenvalue problem which can be solved by 
linear algebraic techniques. The Lie algebraic approach to quasi-exactly solvable problems 
requires that the subspace in question be invariant under a Lie algebra of differential 
operators g, in which case the Hamiltonian belongs to the universal enveloping algebra of q; 
see [16, 17, 37, 47, 52] for details and [23, 29] for applications to molecular spectroscopy, 
nuclear physics, and so on. Zhdanov, [56], indicated how one could characterize quasi- 
exactly solvable operators using higher order symmetry methods. We should also mention 
Hel-Or and Teo, [21], who have applied group-invariant subspaces in computer vision, 
naming their elements "steerable functions" . 

Motivated by a problem of Bochner, [4] , to characterize differential operators having 
orthogonal polynomial solutions, Turbiner, [44] , initiated the study of differential operators 
leaving a polynomial subspace invariant. In one-dimension, the remarkable result is that 
the operators leaving the entire subspace of degree < n polynomials invariant are the 
quasi-exactly solvable operators constructed by Lie algebra methods. These results were 
further developed for multidimensional and matrix differential operators, and difference 
operators by Turbiner, Post and van den Hijligenberg, [39,40,45,46], and Finkel and 
Kamran, [8]. 

In this paper, we broaden the general theory of nonlinear separation to include partial 
differential operators, and argue that it constitutes the proper nonlinear generalization of 
quasi-exactly solvable linear operators. Our theory provides an explicit characterization of 
all nonlinear differential operators that leave a given subspace of functions invariant. In the 
time-independent case, solutions lying in the subspace are obtained by solving a system of 
nonlinear algebraic equations. For evolution equations and certain other dynamical partial 
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differential equations, we show how exphcit solutions in the given subspace are obtained by 
reducing the partial differential equations to a finite-dimensional dynamical system. We 
illustrate our method with a number of significant examples. 

Our methods can also be compared and contrasted with the more algebraic theory of 
D-modules, cf. [2, 7, 31, 32]. In particular, we describe new algorithms for determining 
the annihilator of a given finite-dimensional subspace, based on the study of generalized 
Wronskian matrices and their ranks. These methods derive their justification from the 
general theory of prolonged group transformations developed in [36, 38]. 

The first section of the paper outlines the basic setting of our methods — finite- 
dimensional spaces of analytic functions defined on an open subset of real Euclidean space. 
The case of analytic functions of several complex variables is more subtle, and requires 
cohomological or geometric restrictions on the domain. Section 3 presents the basic tools 
in our study: a multi-dimensional generalization of the classical Wronskian condition for 
linearly independence of functions. In fact, these results form a very particular case of a 
general study of orbit dimensions of prolonged group actions formulated in [38]. The key 
concept is the notion of a "regular" space of analytic functions, and only for such sub- 
spaces is one able to characterize the differential operators that annihilate the subspace 
or leave it invariant. While the Hilbert basis theorem is not generally applicable to ideals 
of analytic differential operators, one can nevertheless algorithmically determine a finite 
generating set of annihilating differential operators when the subspace is regular. Section 
4 contains our basic approach to the annihilator, which relies on a useful extension which 
we name an "affine annihilator" , which is a differential operator that maps every function 
in the subspace to a constant. The construction ultimately rests on an interesting lemma 
characterizing analytic solutions to a system of variable coefficient linear equations of con- 
stant rank; it is at this critical point that the distinction between the real and complex 
analytic situations becomes evident. We include several examples illustrating our con- 
struction, including a convenient characterization of the affine annihilators of simplicial 
subspace of monomials. Section 5 applies these constructions to characterize all linear and 
nonlinear differential operators that leave a given regular subspace invariant. In the scalar 
case, Svirshchevskii, [49, 50, 51], characterized these differential operators using general- 
ized symmetries, and we prove that Svirshchevskii's symmetry operators coincide with our 
affine annihilators, thereby establishing the generalization of Svirshchevskii's methods to 
analytic functions of several variables. Finally, section 6 outlines how our results can be 
applied to the construction of explicit solutions to linear and nonlinear partial differential 
equations based on the method of nonlinear separation. 

2. Function Spaces. 

Let X C M"^ be an open, connected subset of Euclidean space, with coordinates 
{x^, . . . , x^). Our basic set of allowable functions will the space JF = C^(X) of analytic 
real-valued functions /: X — > M. We may regard JF, depending on the circumstances, as 
either a real vector space, or as an algebra over the reals. 

Even though we shall primarily restrict our attention to real domains and real analytic 
functions, much of the exposition can be adapted to other situations. For instance, a much 
easier situation is when = is the space of all germs of analytic functions at a 
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single point a = {a^, . . a"^) e X, or, equivalently, the space Clx^ — a^, . . — a*"]] of 
convergent power scries at the point a. Another interesting example, studied extensively 
in D-module theory, is when JF is the field of meroniorphic functions on X, or, more 
generally, any differential field, cf. [27]. The case when = consists of smooth 

functions on X is also quite interesting, but much more diflBcult to treat owing to a number 
of pathologies that do not appear in the analytic context. With the proper restrictions, 
our methods and results can also be made to apply to various function spaces in the 
complex-analytic category. We briefly indicate how this may be done in Section 4, below. 

Once we have fixed the function space our primary object of study are finite- 
dimensional subspaces M. d In particular, given functions fi,...,fj^ e J^, we let 
A4 = {/i, • • ■ , /fc} denote the subspace spanned thereby. We shall use r to denote the 
dimension of M., and often use f^,...,fj.&Tto denote a basis. In applications to quasi- 
exactly solvable quantum problems, is a finite dimensional module for the action of a 
transformation group G on X, cf. [16, 17]. 

Let V = P(jF) denote the space of linear differential operators whose coefficients 
belong to the function space JF. The multiplication of differential operators is defined in 
the usual manner, making V into an associative, but noncommutative algebra over — 
the latter acting by left multiplication. Furthermore, the operators in V define linear maps 
L: T ^ T , and so T will also be regarded as a P-module. There is a natural filtration on 
D, with "D" denoting the jF-submodule consisting of differential operators of order < n. 
The domains X considered here are such that every linear differential operator L e T>^ is 
realized as a finite sum 

L = Y,hj{x)dj, where dj = , (2.1) 

and the coefficients hj G JF. The sum in (2.1) is over symmetric multi-indices / = 
(ii,...,ifc) of orders < /c = #/ < n. Since = ■ X>", we may identify the 

quotient /V'^~^ with the .F-module of homogeneous differential operators of order n. 

We note that there are = ("^^"~"^) different symmetric multi-indices / of order 
#/ = n, which is the number of different n^^ order partial derivatives dj. Similarly, there 
are 

,(n)=,„ + ,^ + ... + ,„=(™ + ") (2.2) 

different symmetric multi-indices / of order ^/ < n. An n^^ order differential operator 
is uniquely determined by its q^"^^ different coefficients hj{x), < n. Therefore, we can 
identify T>^ with the space JF^"?' ' — the Cartesian product of q^"^^ copies of JF. Explicitly, 
the isomorphism a^-.V^ — > maps a linear differential operator (2.1) to the column 

vector 

cTjL) = h(x) = (...,/i,(x),...)^ (2.3) 
whose entries, indexed by the multi-indices of order #/ < n, are the coefficients of L. 
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3. Wronskians and Stabilization. 



The most basic necessary and sufficient condition for the hnear independence of solu- 
tions to a homogeneous hnear scalar ordinary differential equation is the nonvanishing of 
their Wronskian determinant, cf. [20]. Many of our results will rely on a significant multi- 
dimensional generalization of this classical Wronskian lemma, which can be applied to any 
collection of analytic functions. We refer to [38] for details on the following definitions 
and results, including extensions to both smooth and vector-valued functions. 

Definition 3.1. The n*^ order Wronskian matrix of the functions J^, . . . , G JF, is 
the r X q^'^^ matrix 



Wjx) = WJA,...,/J(a;) 



/A(x) ... d,U{x) ... \ 



\ax) ... d,f,{x) ... I 



(3.1) 



whose entries consist of the partial derivatives of the /^'s with respect to the x"s of all 
orders < #/ < n. 

In the scalar case X C M, the standard Wronskian determinant coincides with the 
determinant of the (r — 1)^* order Wronskian matrix W^_i = W^_i [/j^, . . . , /j.], which 
happens to be a square matrix. 

To this end, we define the Wronskian matrix rank function 

pJx) = rankWJA,...,/J(a;). (3.2) 

Note that p^i^x) only depends only on the subspace M. — {/i, . . . , /^} spanned by the given 
functions, and not on the particular generators or basis. In particular, < p^ix) < dim Al. 
Moreover, p^{x) is lower semi-continuous: if Pni^o) = then p^{x) > k for all a; in a 
sufficiently small neighborhood of Xq. The generic Wronskian rank of order n is 

p* = max{p^(a;) | x e X} . 

The first key result is the following: 

Theorem 3.2. Let M. <Z T he an r -dimensional subspace and p*, n > 0, be the 
sequence of its generic Wronskian ranks. Then p* = r = dimAl for n » sufficiently- 
large. Moreover, if wc define the (generic) stabilization order s = min{n | p* = r} of Al to 
be the minimal such n, then we find 

p* < p* < • • • < p*_i < p* = p*+i = p*+2 = . . . = r = dimM. (3.3) 

In other words, once the generic Wronskian ranks become equal, they stabilize, and 
achieve a value equal to the dimension of the subspace. In particular, the ranks cannot 
"pseudostabilize" , [36, 38], and so must strictly increase before stabilization sets in. It is 
entirely possible, though, that they increase only by 1 at each order — an evident example 
occurs when all the functions only depend on a single variable. 

While knowledge of the stabilization order reduces the amount of work required to 
compute the order of the generic Wronskian rank, one can, in all cases, replace s by r — 1. 
Thus, one has: 
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Corollary 3.3. The analytic functions f-^^{x),...,f^{x) are linearly independent 
if and only if the generic rank of their (r — 1)®* order Wronskian matrix W^_j^(a;) = 
W^-i [/d • • • ) fr] (^) ^'s equal to r. 

Proof: Indeed, Theorem 3.2 implies that /j^,...,/^ are linearly independent if and 
only if p* = r for any n > s greater than or equal to the stabilization order. Moreover, 
(3.3) implies that the stabilization order of an r-dimensional subspace Ai G T is always 
bounded by s < r — 1. Thus, one typically does not need to compute the rank of the order 
n = r — 1 Wronskian to detect linear independence — checking at the stabilization order 
is suflBcient. Another way of stating this result is the following: if the generic ranks of the 
Wronskians W^(x) and Wj^_,_j^(a;) are equal, then this rank is the same as the dimension 
of A4; moreover, this equality first occurs when k = s < r — 1. Q.E.D. 

Kolchin, [27; p. 86], states and proves a version of Corollary 3.3 assuming that the 
functions belong to a differential field, e.g., the field of meromorphic functions. It would 
be of interest to adapt our constructions to this case. 

It is noteworthy that Theorem 3.2 is a particular case of a general theorem governing 
the orbit geometry of prolonged transformation groups acting on jet bundles, described in 
detail in [36]. In the present situation, the relevant group is the elementary r-parameter 
abelian group 



acting on i? = X X M. The dimension of the prolonged group orbits contained in the n^^ 
order jet fiber ^^E\^ equals the Wronskian rank p^{x). Moreover, in the analytic category, 
the group (3.4) acts effectively if and only if the subspace A4{fi, ■ ■ ■ , f^} has dimension 
r. The extensions of this result to more general smooth group actions are also treated in 
[38], and can be applied to the more subtle case of subspaces A4 C C°° containing smooth 
functions. For brevity, we shall refrain from discussing this more complicated case here. 

Although, for an r-dimensional subspace Ai, the (r — 1)®^ order Wronskian has generic 
rank r, it may certainly have lower rank at particular points in the domain X. A simple 

2 

example is provided by the functions fi{x) — 1, f2ix) = , which have first order 

2 

Wronskian determinant detWj^(x) = 2xe^ , which is singular at a; = 0. (The classical 
Wronskian lemma implies that these functions cannot be common solutions to a regular, 
homogeneous, second order linear differential equation.) Our applications typically require 
that, for some n sufficiently large, the Wronskian rank p^{x) be equal to r at every x e X, 
and so we need to determine when this occurs. In the preceding example. 



has rank 2 everywhere, as do all the higher order Wronskian matrices. Therefore, it is of 
interest in understanding, not just how the generic Wronskian ranks behave, but also how 
the Pn{x) behave at a single point. 




(3.4) 






2xe 
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Theorem 3.4. If dimA4 = r and x G X is any point, then there exists a finite 
integer n = n{x) such that Pj^{x) = r for aU k >n. We call n the Wronskian order of the 
point X, and denote it by 

c<;(a;) = min{n I p^{x) — r} . (3-5) 

fc 

For example, if we take fi{x) = 1, f2{x) = , for A; > 2 an integer, then the 
Wronskian order of a; = is equal to a;(0) = k. Therefore, the Wronskian order of a point 
can be arbitrarily large. 

Again, Theorem 3.4 is a corollary of a more general theorem about the geometry of 
prolonged transformation groups. In the language of [38], a point at which p^(a;) < r for 
all k would be known as a "totally singular point" for the associated transformation group 
(3.4). However, Theorem 6.4 of [38] states that an analytic transformation group admits 
no totally singular points. In fact, for the elementary group (3.4), this result is not hard 
to prove directly — it is a consequence of the basic result that a nonzero analytic function 
cannot have all zero derivatives at a single point. And this is the main reason why we 
must restrict out attention to analytic, as opposed to smooth, functions. 

Given an r-dimensional subspace At = {fi, . . . , /^}, Corollary 3.3 implies that most 
— meaning those belonging to a dense open subset of X — points have Wronskian order 
r — 1 or less.^ The exceptions are the points x where W^_j^(x) has less than maximal 
rank. This permits us to formulate a basic estimate on the Wronskian order of the space 
A4. In general, given an ordered r-tuple I = {I^, . . . , I^) consisting of symmetric multi- 
indices, we define its order to be = max{^/j^ | = l...r}. Let Wj(x) denote the 
associated r x r submatrix of yVj^{x) whose columns are indexed by the multi-indices in 
I. More prosaically, the entries of Wj(x) are the partial derivatives dj f^, |U, z/ = 1, . . . , r, 
indicated by I. Note that Wj(a;) depends on at most #1*^ order derivatives of the functions 
f-^,...,fj,. The Wronskian minor is then 

Mi(a;) = detWi(a;). (3.6) 

Clearly, the Wronskian order of a point x is the smallest multi-index order with nonvan- 
ishing Wronskian minor, i.e., 

a;(a;) = min{#I I Mi(a;)^0}. (3.7) 

Given a non-zero analytic function h{x) ^ 0, we define its vanishing order 'ccrf^{x) 
at a point to be the order of its first nonzero terms in the power series expansion at x. 
Alternatively, we can define the vanishing order by differentiation: 

Wf^(x) = min{#/ I djh(x) ^ 0} . (3.8) 

In particular, tx7^(x) = if h{x) ^ 0. Note that the vanishing order of a nonzero analytic 
function is always finite. 



^ If dim X = 1, then the Wronskian order is never less than r — 1. More generally, if dim X = m, 
then the minimal Wronskian order I of any subspace is bounded from below by the inequality 

^(0 = (jn+l^ > r. 



7 



Definition 3.5. Given a subspace A4, let Wi{x) = ^^^^^{x) denote the vanishing 
order of its 1*^ Wronskian minor Mj(x). The Wronskian vanishing order of A4 is defined 
as the minimal vanishing order of the minors of 'W^_i{x), so 

w{x) = min {w^ix) \ #1 < r - 1} . (3.9) 

The Wronskian vanishing order of a subspace provides an immediate bound on its 
Wronskian order at a point. 

Proposition 3.6. Let dim = r and let x E X. If the Wronskian vanishing order 
of M. at X is w{x), then the Wronskian order of x is bounded by 

uj{x) <r -l + tu{x). (3.10) 

Proof: It suffices to note that if J is any multi-index of order = j, then the J*^ 
order derivative of a Wronskian minor Mj(a;) of order = i can be written as a finite 
linear combination of Wronskian minors 

djm^{x) = Y, MK(a;) 

K 

of orders #K < i + j. Hence, if J and I are such that 5jMj(a;) ^ 0, then there must exist 
a K with #K <i + j such that M.j^{x) 0. Q.E.D. 

Example 3.7. Consider the three-dimensional subspace Al = {1, cos cos 2a;}. The 
Wronskian determinant of these three functions is det W2(a;) = — 4(sina;)^, and hence M. 
has Wronskian vanishing order 3 at x = titt, n G Z. The fourth order Wronskian matrix 
^^{x) is found to have rank 3 for all a; e M, and hence the Wronskian order at the singular 
points X — nn is, 4. 

If Xq has Wronskian order n, then, by continuity, p^{x) = r for all x in some neighbor- 
hood of Xq. However, a global bound on the Wronskian order may not exist. We therefore 
introduce the following important definition. 

Definition 3.8. A subspace Ai is called regular if it has uniformly bounded Wron- 
skian order at each point, i.e., there exists a finite n such that p^{x) = r for all x E X. 
The minimal such n is called the order of Ai. 

Many of the basic results in this paper require the underlying regularity of the sub- 
space. 

Example 3.9. Not every subspace is regular. For example, let X = M. Consider the 
one-dimensional subspace spanned by an analytic function f{x) which has a zero of order 
k at Xf^, for /c = 1, 2, 3, . . . , where ^ oo as /c — » oo; such a function can be constructed 
using a Weierstrass product expansion, cf. [1; p. 194]. Then rankW^[/] (a;^) = for n < /c, 
while rank [/] (a;) = 1 for n > A; and x in any neighborhood of Xj^ that does not contain 
^k+n^k+2j ■ • • ■ other words, the Wronskian order of each Xj^ is equal to k. 
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We know two straightforward mechanisms for proving that a given subspace is regular. 

The key is to avoid an infinite sequence of points whose common Wronskian order is 
unbounded, as in the preceding example. The first method relies Proposition 3.6 for 
proving regularity. 

Proposition 3.10. IfAi is a subspace with uniformly bounded Wronskian vanishing 
order, so zu{x) < n for aU x E X , then X is regular. 

Alternatively, one can generate regular subspaces by restricting the domain of the 
functions and using the basic properties of compactness. 

Proposition 3.11. If M. C C'^(X) is any subspace, and Y X is an open subset 
with compact closure in X , then the subspace Xi = A4 lyC. C^(F) obtained by restriction 
to Y is a regular subspace. 

4. Annihilators. 

Our characterization of differential operators preserving spaces of analytic functions 
relies on several intermediate constructions of independent interest. The first order of 
business is to characterize the differential operators that annihilate all the elements of our 
subspace Ai. We begin with the linear annihilating operators; their nonlinear counterparts 
will be treated in the following section. 

Definition 4.1. The annihilator A = A{Ai) of a subspace C is the set of all 
linear differential operators that annihilate every function in Ai, so 

A{M) = {KeV\ K[f] = for aU f e M} . 

We will refer to the elements of A{J^) as annihilating operators. We note that A gV 
is, in fact, a left ideal, since if K[f] = and L e D is any linear differential operator, then 
clearly L[K[f]] = and so L ■ K e A. 

Example 4.2. Let X C R. Let M C C'^{R) be an r-dimensional subspace. A 
classical construction, cf. [54], produces an r^^ order annihilating operator 

= h^{x)d'' + + ■ ■ ■ + h^{x). (4.1) 

Indeed, introducing a basis /i, • • • , of /A, then the conditions K^[fj^] — 0, u = 1, . . . ,r, 
forms a homogeneous system of r linear equations for the r + 1 coefficients of K^. Cramer's 
rule produces the solution 

h{x) = {-ir-'^M,{x), (4.2) 

where M^(x) = det Wq^^ ^{x) denotes the r x r Wronskian minor obtained by 

deleting the k^^ column of the r x (r + 1) Wronskian matrix W^(x). In particular, the 
leading term h^{x) = M^(x) = det W^_^(x) is the classical Wronskian determinant, and 
hence is a nonsingular differential operator if and only if the classical Wronskian never 
vanishes, which implies that M. is regular of the minimal possible order r — 1. 

If h^{x) 7^ is never zero, then the annihilator A is generated by K^. Indeed, to prove 
that every other annihilating operator has the form K = L ■ K^, we note that every linear 
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ordinary differential operator can be written (uniquely) in the form T = Q ■ + R, where 
Q and R G T>^~^ has order at most r — 1. (If T has order < r — 1, then Q = 0.) But 
T G A{Ai) if and only if R E A{Ai) is an annihilating operator. But the dimension of the 
kernel of a linear differential operator of order k is at most^ k, and so i? = 0. 

On the other hand, if h^{x) vanishes at points x & X, then does not provide a basis 
for the annihilator of M.. For example, the annihilator of the one-dimensional subspace 
spanned by the function fi{x) = x is generated by the two differential operators xd — 1 
and 9^. (Although d'^ = ^d- (xd—l), the operator is not analytic, and so not allowed.) 
More complicated cases are discussed below. 

Example 4.3. Let X C M"^ and consider a one-dimensional subspace Ai = {/} 
where / ^ 0. The linear operators — fd^ — /j, i = 1, . . . , m, where = d^f — df /dx^, 
clearly belong to A^. If f{x) ^ never vanishes, then K^^, . . . , form a basis for A. The 
proof of this fact is similar to the ordinary differential operator result of Example 4.2. We 
first note that any differential operator can be written in the form T = ^Q^- K^ + g, where 
Qi, • • • , Qjri ^ ^ differential operators, and g G ^ is a multiplication operator. 
Clearly T G A{A4) if and only if g = 0, proving the result. The case when / vanishes on a 
subvariety of X is more subtle — see below. Extensions to higher dimensional subspaces 
are also discussed below. 

The filtration of T> induces a filtration of the annihilator, and we let A^ = Ar\ 
denote the subspace of annihilating operators of order at most n. 

Proposition 4.4. Let the analytic functions /i,...,/^ span a Gnite-dimensional 
subspace M = {fi, . . . , f^} C J^. Let the associated n^^ order Wronskian matrix. The 
map defined by (2.3) defines an isomorphism 



Proof: It is sufficient to note that for every K G A'^{Ai) and corresponding vector 
h = (7^{K) with analytic entries, the entry of the matrix product • h is equal to 



Had we taken T to be the algebra of power series (or analytic germs), then a straight- 
forward adaptation of the classical Hilbert basis theorem, [22, 55] would prove that the 
annihilator A{M.) of every finite-dimensional subspace M. <Z Tis finitely generated. This is 
because the power scries algebra is Noetherian, cf. [55]. However, the algebra T — C'^{X) 
of globally defined analytic functions is not Noetherian, and so the Hilbert basis theorem 
does not apply. The following example shows that not every ideal of T — C^^X) is finitely 
generated. 



' Operators with degenerate symbols may not admit enough analytic solutions to span a full k- 
dimensional kernel. For example, the functions annihilated by the first order operator L = xdx + 1 
are multiples of the non-analytic function and so L has a zero dimensional (analytic) kernel. 




(4.3) 



between the n^^ order annihilator of Ai and the common kernel 




Q.E.D. 
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Example 4.5. Consider again the analytic function f{x) introduced in Example 3.9. 
Let f^'^\ /c = 0, 1, 2, . . . denote the corresponding derivatives and note that f^^^d — /(^+^) 
is an annihilating operator for all k. It isn't hard to see that the ideal of JF = C^(X) 
that is generated by all the f^^^ is not finitely generated, and therefore, in this case, the 
annihilator A oi M = {/} is not a finitely generated ideal of V. 

Later we will prove that the annihilator ideal of a regular subspace is finitely generated, 
to this end, we introduce a useful extension of the notion of the annihilator. 

Definition 4.6. The affine annihilator K,{}A) d T) oi a, subspace M. d is the 
subspace of those linear differential operators that map every function in to a constant 
function: 

1C{M) = {LeV\ L[f] = c e M, for all / G M} . 

Note that, as defined, the affine annihilator 1C{M.) is a vector space, rather than an 
.F- module. Considered as a vector space, the annihilator A{M.) is evidently a subspace of 
K,{M.). The difference between the two subspaces has a natural interpretation. 

Definition 4.7. The operator dual to the subspace M. is the quotient vector space 

M* = 1C{M)/A{M). (4.4) 

Since we are quotienting by the annihilator, there is a natural action oi M* on Al 
itself induced by the action of /C(A1). In this fashion, if M. is finite-dimensional, then 
there is a natural linear injection from M.* into the abstract dual oi M.. As the following 
theorem will demonstrate, regularity implies that this injection is, in fact, an isomorphism. 

Theorem 4.8. If Ad d T is a regular r -dimensional subspace of order s, then its 
operator dual Ai* is also r -dimensional. Moreover, if f^, . . . , forms a basis of M., then 
there exists a dual basis for M.* represented by differential operators L^, . . . , of order 
at most s such that 

Liifj) = 5}, i,j = l,...,r. (4.5) 
Example 4.9. Let us consider the "triangular" or simplicial polynomial subspaces 

T^ = {xY I 0<i + j<n}. (4.6) 

We remark that forms a module for the standard representation of the Lie algebra sl(3) 
by first-order differential operators, cf. [16], and so plays an important role in the theory 
of quasi-cxactly solvability and orthogonal polynomials, cf. [44, 45] . It is not hard to see 
that 7^ forms a regular subspace of order s = n. 

To construct the dual basis, we look for a set of differential operators L^j, < i+j < n, 
such that 

L..{x''y^) = S^^Sji, foraU < i + j < n, < k + I < n. (4.7) 

Let 

S = xd^ + ydy 
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denote the degree or scaling operator. Let 

Pk{^) = ^-^ {x - l){x - 2) ■ ■ - {x - k) 

be the unique polynomial of degree k such that P;j(0) = 1 and p{j) = for j = 1, . . . , A;. 
Then dual basis operators are given by 

L^, = ^^Pn-i-,{S).^l^l. (4.8) 

Indeed, the two last factors annihilate any monomial of degree < i -\- j except for x'^y^ , 
whereas the polynomial in S will annihilate the higher degree monomials. For n = 2, the 
dual basis is explicitly given by 

-^20 = \^l-> ^11 = ^xy^ -^02 = 

Lio = i-S + 1) • 5. = i-^d, - yd^ + 1) • d,, 

L^^ = (-S+l).dy = i-xd,-ydy + l)-dy, ^ ■ ^ 

L,,= lS'-IS + l^ \x^dl + xyd,^ + ly'dl - xd, - yd^ + 1. 

Formula (4.8) readily generalizes to the simplicial modules in m variables, i.e., the poly- 
nomial subspaces generated by the monomials 

x^ = xY ■■■xl\ ^I = J2ii< n. (4 JO) 

I 

We note that this subspace forms a finite-dimensional module for the standard represen- 
tation of sl{n) by first-order differential operators. 

The proof of Theorem 4.8 ultimately rests on the following technical Lemma. 

Lemma 4.10. Let A{x) be a real-analytic, r x n matrix-valued function, deGned 
for X E X G M"*. Suppose rankA(x) = r for all x E X . Let b:X W be an analytic 
vector-valued function. Then there exists a solution /i: X — > R" to the matrix equation 
A{x)h{x) = b{x) which is analytic for all x. 

Proof: Clearly, it suffices to prove the result when b{x) = is constant, equal to the 
j^^ basis vector. Let hj (x) denote the sought-after solution with this particular choice of 
the right hand side. Let v-^^{x), . . . , v^{x) denote the rows of A, considered as vector- valued 
functions Vf^: X Let *: A M'' ^ A'"" M"" denote the flat space Hodge star operation 

in the exterior algebra A*]R", cf. [53; p. 79]. In particular, *1 = A - ■ ■ Ae„ G A'^M"^ is the 
volume form. One can write the matrix system A{x)hj{x) = e- in the equivalent exterior 
form 

v^{x)A[^h^{x)]^{ °' (4.11) 

obtained by applying * to each equation and using the fact that v ■ w = * {v A *w). An 
evident solution to the first r — 1 equations in (4.11) is 

hj{x)= * [v^{x) A ■ ■ ■ AVj_-i^{x) AVj^-i^{x) A ■ ■ ■ Av^{x) A^{x)] , (4.12) 
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where 7:^ — > A" is an arbitrary analytic map. The final equation in (4.11) leads to 

(-1)"+^' v^{x)A---A v^{x) A 7(x) = * 1. (4.13) 

This is a linear equation in the coefficients gi{x), . . . ,gf.{x) of j{x), where k = (^), and as 
such has the form 

a^{x)g^{x) H + a,^{x)g,^{x) = 1, (4.14). 

Here a]^(x), . . . , aj^{x) are analytic functions explicitly determined by Vi{x)j . . . , t'^(x); in- 
deed, up to sign, they are just the rank r minors of A{x). Moreover, these minors cannot 
simultaneously vanish since Vi{x), . . .,Vj,{x) are assumed to be everywhere linearly inde- 
pendent. Since we are dealing with real-valued functions, an evident solution to (4.14) 
is 



Substituting (4.15) into (4.12) and using the fact that *{luA *a;) = || a; ||^ for any uj e A*I 
produces an explicit analytic solution to (4.11) in the form 

V;^^{x) A • • ■ AVj_;^{x) A Vj^-^^{x) A ■ ■ ■ A v^{x) A*{v^{x) A ■ ■ ■ AVj.{x)) 



) n 



hjx) = {-1)''+^ * 



v-i^{x) A ■ ■ ■ A v^{x) II 



(4.16) 

Note that the denominator is nowhere zero since the rank of A{x) equals r. Q.E.D. 

Lemma 4.10 plays a critical role in our theory. It relies on the fact that there exists 
an analytic solution a-^{x), . . . , af,{x) of (4.14), which, in the real category considered here, 
is constructed in (4.15). The advantage of choosing the reals as the ground field lies in the 
fact that our results are valid for arbitrary domains. However, it is important to keep in 
mind that the solutions to equation (4.14) given by (4.15) may fail to be analytic in the 
complex category, and indeed. Lemma 4.10 is not true for all complex domains X. 

Example 4.11. A simple counterexample is provided by the case X = C'^ \ {(0, 0)}, 
and the 1x2 matrix function A{x,y) — {x,y), which has rank 1 for all {x,y) G X. There 
is no complex-analytic vector-valued function h{x^y) = {hi{x,y), h2{x,y))'^ such that 

A{x,y) h{x,y) = xh-^{x,y) + yh^ix.y) = 1, (4.17) 

for all {x, y) 7^ (0, 0). This stems from the fact that X is not a domain of holomorphy, and 

so any complex-analytic function on X can be extended to a complex analytic function on 
all of C^, cf. [28; §0.3.1]. Extending h^{x,y) and h2{x,y) in this manner, by continuity 
(4.17) would also have to hold aX x — y — 0, which is clearly impossible. 

The preceding example makes clear that the complex case is more subtle, and requires 
additional assumptions. One could, for instance, restrict oneself to the case where JF is the 
algebra of complex- analytic germs (i.e., convergent power series). In this case, the natural 
analogue of the order of an r-dimensional module would be the smallest n such that the 
rxn matrix formed by the constant terms of the n*^ order Wronskian W has rank r. For 
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such an n, one of the rank r minors of would be a unit, thereby making Lemma 4.10 
true. 

The utiUzation of more general domains X C C"^ would require a cohomological 
assumption in order that (4.14) admit an analytic solution. One could, for instance, 
demand the vanishing of every whose coefficients lie in a coherent sheaf, or more 
generally that X be a Stein manifold, [19]. Alternatively, one could impose a geometric 
restriction on the domain X. For example. Corollary 7.2.6 in Krantz, [28; p. 293], states 
that if X C C"^ is a pseudoconvex subdomain, then (4.14) has a complex analytic solution 
Qiix), . . . ,gk{x) provided the a^^'s do not simultaneously vanish. Therefore, Lemma 4.10, 
and its consequences, would hold provided X is pseudoconvex. To keep matters simple, 
though, we shall not return to the complex- analytic situation. 

Proof of Theorem 4-8: Recall first that we are seeking dual basis operators of or- 
der s. To construct an i^^ dual basis operator, we use the isomorphism (2.3) to rewrite 
equation (4.5) in the equivalent matrix form 

W,(x)h,(x)=e,, (4.18) 

where h.-{x) = crg{L-) E T^'^ is the vector of coefficients of the operator L^, and is 
the standard i*^ basis vector for W . Since s is the order of the subspace M.^ we have 
rankW^(x) = r for all x^ and hence Lemma 4.10 shows that there exists a solution h^(a;) 
to (4.18) which is analytic for all x E X. This simple construction produces the required 
basis for the affine annihilating operators. Q.E.D. 

Having proven the existence of a dual operator basis, we can return to the examination 
of the annihilator ideal. Thus, for the remainder of this section we assume that is a 
regular r-dimensional subspace of functions with basis fi,...,f^, and fix a dual basis 
L^, . . . G }C{A4) of operators of order s or less. It is important to keep in mind that 
M*, as defined, is a vector space, and not an J^-module. When the need arises, we will use 
Al* to denote the .F-module generated by L^, . . .,L^. Our results rely on the following 
key observation, whose straightforward proof is left to the reader. 

Proposition 4.12. The mapping a:T> ^ A given by 

r 

a{T) = T-Y,T[f.]L,, TeV, (4.19) 

defines an T-module homomorphism. Furthermore, a is a left inverse of the inclusion 
homomorphism i: A^T>, thereby yielding the -module decomposition D ~ ^ eAi*. 

At this point it is also important to note that a basis f-^^,...f^ofAi does not uniquely 
determine the dual basis L^, . . . , L^; one can obtain other bases by adding elements of A'^ 
to the operators L^. For this reason the projection a:V ^ A is not a natural object, but 
rather depends on the choice of the dual basis. 

Remark: On the other hand, for a simplicial module 7^ described in Example 4.9, 
there are no annihilators of degree less than n + 1 = s + 1, and hence the dual module 
Ai* is uniquely determined, or equivalently, A^ = 0. Thus, the splitting of D ~ ^ ®M.* 
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described in Proposition 4.12 is canonical in this case. Subspaces having this property will 
be called saturated, and form an interesting class worth further investigation. 

For each multi-index / let us denote the "basic" annihilator 

Kj = a{dj) (4.20) 

obtained by projecting the basis differential operators dj onto the annihilator A. It is 
important to note that since the have orders at most s, if #/ > s + 1, then the leading 
order term (or symbol) of Kj is just dj. 

Corollary 4.13. Any differential operator T & T> can be written as a finite linear 
combination 

r 

T = Y1 9M K + Y. ^/(^) (4-21) 

of the dual affine annihilators and the basic annihilators. 

Proof: Since the operators dj^ generate V as an jF-module, we may use Proposi- 
tion 4.12 to infer that L^, . . . along with the operators Kj also generate V ; and that 
the affine annihilators along with all Kj such that #/ < s generate T>^. It follows imme- 
diately that for n > s, the Kj such that #/ < n generate A^. Q.E.D. 

Remark: The linear combination (4.21) is not necessarily unique. However, one can 
eliminate precisely r of the Kj's in order to suppress any redundancy and thereby yield a 
well-defined .F-module basis {L^, Kj} for T>. 

Theorem 4.14. Let M. <Z T be a regular r -dimensional subspace of order s. Then 
its annihilator ideal A{A4) is finitely generated by differential operators of order at most 
s + 1. 

Proof: We shall prove that the operators Kj such that #7 < s -|- 1 generate A as an 
ideal of V. Clearly all of A^~^^ can be so generated. Let J be a multi-index whose order 
is > s -|- 1, and choose multi-indices /, with ^/ = s + 1 and such that dj = d^dj. 
According to the remark at the end of the proof of Corollary 4.13, K j and dj^ ■ Kj have 
the same leading term, and therefore the order of the difference Kj — dj^ ■ Kj of the two 
operators will be smaller than the order of J. The desired conclusion now follows by 
induction. Q.E.D. 

Remark: The set of generators Kj constructed above is typically not minimal. How- 
ever, since the set in question is finite, minimal generating sets do exist. 

Remark: In the theory of D-modules, one is interested in subspaces generated by 
rational functions f{x) = p{x)/q{x), with polynomial p,q. However, the annihilating 
operators of interest are required to have polynomial coefficients, and so the set-up is a 
bit different from that considered in this paper. See [32] for applications of powerful 
techniques from Grobner bases and the theory of -D-modules towards the determination 
of the annihilators of such rational subspaces. It would be interesting to see whether our 
techniques have anything to add to this theory. 
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Example 4.15. Consider the two-dimensional subspace A4 = {l,cosa;}. The second 
order Wronskian matrix is 



1 

cos X — sin X — cos x 



and hence A4 is regular of order 2. Applying the algorithm of Lemma 4.10 leads to the 
dual operators 

= cos^ xdl + cos a; sin a; + 1, = — cos xd^ — sin x d^. 

It's important to mention that often the algorithm does not produce the most efficient 
answer. Indeed, for this particular M., a more stream-lined basis of dual operators is given 

by 

L-^^ = ^'^ + l, L2 = — cos X 9^ — sinx 9^. 
The annihilator ideal, A{Ai) is generated as a ring by the following operators: 

^1 = -sinxS^ + cosxSg., A2 = dl + d^. 

This is a consequence of Theorem 4.14, which states that A{M) is generated by the 
projections of the order 3 basic differential operators via (4.19). Indeed the projected 
operators are: 

a(l) = sinx a{dy.) = cos X A-^^, a(9^) = — sinx ^i, a{dl) = A2 — cos x A-^^ 

Therefore A-^ and A2 suflBce to generate all of A{M.). 

Example 4.16. The subspace M. = {1, cos cos 2a;} considered in Example 3.7 is 
considerably more difficult, since one needs to use the fourth order Wronskian to implement 
the algorithm. The explicit formulae are quite complicated. Set 

iy{x) = 144 - 432 sin^ x + 432 sin^ x + 400 sin^ x. 

Indeed, ^{x) is the sum of the squares of the rank 3 minors of the fourth order Wronskian. 
The reciprocal of ^{x) will therefore be a factor in the expressions for the dual operator 
basis obtained from the algorithm of Lemma 4.10. 

The dual operators are given by 

i'(x)L-^^ = i/(x) + (180 cos a; sin a; — 612 cos x sin^ x + 696 cos x sin^ x)dj. + 
+ (180 - 720 sin^ a; + 1164 sin^ x - 616 sin^ x)dl + 
+ (36 cos a; sin a; + 180 cos x sin^ x — 456 cos x sin^ a;)5| + 
+ (36 - 144 sin^ a; - 60 sin^ x + 136 sin^ x)d^, 

v{x)L2 = (-192 sin x + 656 sin^ x - 736 sin^ x)d^ + 

+ (—192 cos X + 576 cos x sit? x — 656 cos x sin^ x)d'^ + 

+ (-48 sin X - 176 sin^ x + 496 sin^ x)dl + 

+ (—48 cos X + 144 cos x sit? x + 176 cos x sin^ x)d^, 

v{x)L^ = (12 cos a; sin a; — 20cosa;sin^a;)53. + (12 — 24sin^a; + 20sin^a;)5^ + 

+ (12 cos a; sin a; + 20 cos x s\t? x)dl + (12 - 24 sin^ a; - 20 sin^ x)dl, 
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Once again, since we are not dealing with a saturated module, the choice of the dual 
operators is not canonical, and there exists a more stream-lined dual basis: 

L2 = — I cos x 9^ — sin a; 9^ — I cos a; 9^ — | sin x 

L3 = + I sin X cos X 9^ + + I sin^ x) 9^ — | sin x cos x 

The annihilator ideal, A{M.) is generated as a ring by the following operators: 

= 4(cos^ X - l)dl + 1 2 sins cos ,2; - 4(2cos^a: + 1) 
A2 = — sin xd'^ + cos xd^ — 4 sin xd'^ + 4 cos x 
A^ = dl + 5dl+4d^ 

To confirm this one needs to check that A^,A2,A^ generate the projections ci(9^), n = 

0. . . . , 5. Indeed the projected operators are given by: 

a(l) = — ^ sin2a;^i + | sina;^2; '^(^x) — — ]^ cos2a;^i, 

0(5^) = I sin 2xA-^ + ^ sin x A^., ^ip^^ = \ cos 1x A-^ + cos x A^., 

0(5^) = — I sin2a;^i — | sina;^2; "^(^a) — — | cos2a;^i — 5cosa;^2 + ^3 

5. Nonlinear Operators and Invariant Subspaces. 

Let us next turn to the characterization of nonlinear annihilating and affine annihi- 
lating operators. We let E = S{X) denote the space of nonlinear differential operators. 
More specifically, in the case of analytic functions, = C^(X), the space £ consists of 
all analytic differential functions, cf. [36], meaning analytic functions F:J"'(X, R) M 
defined (globally) on the n^^ order jet space of real-valued functions on X. Here n > is 
the order of the differential function F. We write F[u] = F{x, w*^"-*) for such an operator, 
where the square brackets indicate that F depends on x,u and derivatives of u. 

Our main goal is to prove a structure theorem for those differential operators, both 
linear and nonlinear, which preserve a given subspace A4 C J^. Throughout this section 
we assume that is a regular r-dimensional subspace of order s, with basis /i, • • • , and 
dual basis L^, . . . , L^. Let ^/{Ai) C S denote the set of nonlinear annihilating operators, 

1. e., those operators that map every function in Ai to zero. Using our basic annihilating op- 
erators (4.20), we can readily construct the most general nonlinear annihilating differential 
operator for our subspace. 

Theorem 5.1. Every operator F e J\f{M.) can be written as a Gnite sum 

F[u]^Yl GM-Kj[u], (5.1) 

where the Gj are arbitrary elements of £ and Kj = a{dj) are the basic annihilating 
operators. 
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Proof: From Corollary 4.13, we know that there is an analytic function G{x,y,z), 
where y = (y^, . . . , y^), z = { . . . , Zj, . . such that 

F[u]=G{x,L,[ul...,LM,---,Kj[u],...), (5.2) 

depending on finitely many of the Kj. If it = /j(a;) G Ad, then L^[u\ = q e M, and so 
substituting u into (5.2) yields 

F[u] = G(x,c,(i) = 0, where c = (c^, . . . , c^). (5.3) 
Equation (5.3) will hold for all x, c if and only if 

Gix, y,z) = Y^ Gi{x, y, z) z^. 
I 

Therefore, 

q 

G{x, LiM, . . . , L,M, . . . , K,[u]i- • • ) = E ^/M' 

where the coefficients Gj[u] are differential functions. Q.E.D. 

Let K^,...,Ki be a (minimal if desired) generating set for the annihilator A. For 
example, according to Theorem 4.14, one can choose the from among the basic anni- 
hilators Kj of orders < s + 1. Then we can immediately simplify (5.1) to only use the 
generating annihilating operators. 

To be precise, let Q = Q{X) denote the ^-module consisting of differential operators 
whose coefficients are differential functions. Such an operator Z e ^ is given by a finite 
sum 

Z = J2Gj[u]d, = J2 Gj{x,u^-^) d, . (5.4) 
/ I 

Note that the operator G:£ ^ £ maps differential functions to differential functions. 
Corollary 5.2. Every nonlinear annihilating operator F e N'{M.) has the form 

I 

F = J2 Z^-K, (5.5) 

where Z-^^, . . . ,Zi are arbitrary elements of Q. 

We are now in a position to prove the main result of the paper. Let 

V{M) = {P&V\ P[M\ c M] , Q{M) = {QeE \ Q[M\ c M] (5.6) 

denote, respectively, the left ideals consisting of all linear, respectively nonlinear, differen- 
tial operators that preserves the given subspace M.. 
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Theorem 5.3. Let M. d T he a regular r -dimensional subspace of analytic func- 
tions of order s. Let L^, . . . , G IC{A4) be a dual basis for its afRne annihilator, and 
K^, . . . , Ki G A{Ai) a generating set of annihilating operators. Then every nonlinear 
operator Q G Q{A4) that preserves M. can be written in the form 

r I 

Q[u]=J2 f^i^)H^{LM,■■■,LM) + Y. ZM-KM, (5.7) 

i=l v=l 

where the G C^(R'^) are arbitrary analytic functions, and where the G Q{X) are 
arbitrary operators. 

Proof: Suppose f{x) = ^Cjfj{x) G M.. Since Q[f] G M., for each z = l,...,r we 
have Lj(5[/]] = H^{c-i^, . . . ,c^) is a constant depending on the coefficients of /. It foUows 
immediately that 

r r 

E E fiix)H,{L,[u],...,LM) = Q 

i,j=l i=l 

for aU u G A4, and so the result (5.7) follows immediately from (5.5). Q.E.D. 

In particular, every linear operator P G V{M.) that preserves M. can be written in 
the form 

r I 
i,j=l u=l 

where the a^^ G R are arbitrary constants, and the R^^ &V are arbitrary linear differential 
operators. Therefore, we conclude: 

Corollary 5.4. The space of linear differential operators leaving a regular subspace 

M. (Z T of analytic functions invariant can be decomposed into a semi-direct product 
V{M) ^ End(A^) K A{M) of the space End(>l) ^ {A: M ^ M] ~ M®M* of linear 
endomorphisms with the annihilator ideal A{M.). 

In the case when M. is generated by polynomials, this corollary is proved directly by 
Post and Turbiner, [39]. Thus, the class of regular subspaces form the proper analytic 
generalization of the polynomial subspaces. See also [8, 40]. 

In the one-dimensional case, the expressions L^[f] are directly related to the gener- 
alized symmetries considered by Svirshchevskii, [49]. The following example serves to 
illustrate this. 

Example 5.5. Let p = 1 and let M. denote the space of quadratic polynomials 
ax^ + 6x + c in the scalar variable x, with basis x^. The annihilator is generated 
hj K = since the subspace forms the solution space to the linear ordinary differential 
equation 

= 0. (5.9) 

According to [49; Example 3], the fundamental invariants for (5.9) are I^{x, u^"^^) = Ji,[u\, 
where 

= dl, J2 = xdl - d^, J3 = x'^dl - 2xd^ + 2. 
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Since 

J,[1]=0, J2[l] = 0, J3[l]=2, 

Jl[x] = 0, -^2M — ^si^] = 

Ji[x2]=2, J2[x2]=0, J3[x^]^0, 

we see that L^^ = | J3, L2 = —J2, -^3 = |<^i ^re the dual operators with respect to the given 
basis oi Ai. In accordance with the general result, every nonlinear differential operator 
Q G Q{A4) leaving A4 invariant takes the form 

Q[u] = Aq[u] + A^[u]x + A^iulx^ + T[u\ ■ K[u], 

where 

with e C'^(M^) arbitrary, and where T e ^ is an arbitrary (nonlinear) differential 
operator. 

The precise connection between our approach and that of Svirshchevskii relies on the 
theory of generalized symmetries of differential equations, as presented, for example, in 
[34; Chapter 5]. We assume that the reader is familiar with this theory for the remainder 
of this section, and take tt to be a scalar variable, although vector-valued generalizations 
are straightforward. 

Proposition 5.6. Let A[tt] = be an analytic, homogeneous system of linear differ- 
ential equations and let M. <Z T denote the vector space of solutions. Then the generalized 
vector field Vg = Q[u] is a symmetry of A = if and only if the differential operator 
Q[u\ e Q{M.) leaves M. invariant. 

Note that Q[u] may be nonlinear, although many nondegenerate linear partial differ- 
ential equations only admit linear generalized symmetries; sec [48]. Proposition 5.6 is an 
immediate consequence of the basic definition of generalized symmetry, which requires that 
the vector field Vq leave the solution space to A = infinitesimally invariant. Linearity of 
the differential equation implies that infinitesimal invariance coincides with invariance of 
the solution space. Thus, if the differential equation is of finite type, [18], its solution space 
is finite-dimensional, and hence we can completely characterize the generalized symmetries 
of the system using the affine annihilators of the solution space along with Theorem 5.3. 

In Svirshchevskii's method, one considers an r-dimensional subspace Ai consisting of 
analytic functions of a single variable x. Assume, for simplicity, that Ai is regular of order 
r — 1. The linear ordinary differential equation that characterizes A4 is given by 

KM] = 0, (5.10) 

where is the r^^ order differential operator (4.1) constructed in Example 4.2. The 
resulting formulae (5.7) for the generalized symmetries Vg of (5.10) is written in [49] 
in terms of the first integrals of (5.10). Recall that a differential function F{x,u^'^'^^^) 
depending on at most (r — 1)^*^ order derivatives of u is called a first integral of the r^^ order 
ordinary differential equation (5.10) if and only if its derivative D^F = vanishes on the 
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solutions. The first integrals of a linear ordinary differential equation can be constructed 
in terms of the solutions to the adjoint equation 



(5.11) 



where K* is the usual (formal) adjoint differential operator, of. [34; p. 328]. Indeed, inte- 
gration by parts shows that if z{x) is any solution to the adjoint equation (5.11), then 



determines a first integral to (5.10). The explicit formula for the first integral F[z,u] is 
classical, and described in [49]. 

On the other hand, since a first integral is constant on the solution space A4 to 
(5.10), it is an affine annihilator. Choosing a basis Zi, . . . , for the solution space to the 
adjoint equation (5.11) produces r linearly independent first integrals, and hence a basis 
for the operator dual of the subspace. Thus, in this manner we recover Svirshchevskii's 
construction of the generalized symmetries of homogeneous, linear ordinary differential 
equations. 

Conversely, given a nonlinear ordinary differential equation, Q[u] = 0, one can con- 
struct the subspaces Ai G T it leaves invariant by the following "inverse symmetry proce- 
dure" . One needs to classify all linear ordinary differential equations which admit Vg as a 
generalized symmetry. For example, it can be shown that if Q^x^u^"^^) is a nonlinear gen- 
eralized symmetry of a nonsingular linear ordinary differential equation (5.11) of order r, 
then, necessarily r < 2n + l. Consequently, any invariant subspace of a nonlinear n^^ order 
differential operator has dimension at most 2n -|- 1. In practice, the determination of the 
equations of a prescribed order possessing a given generalized symmetry is a straightfor- 
ward adaptation of the usual infinitesimal computational algorithm for symmetry groups 
of differential equations, cf. [34]. Examples of this procedure appear in [49, 50, 51]. 

For functions depending on more than one independent variable, a similar procedure 
works, although now one can no longer use a single linear partial differential equation to 
characterize the subspace, but must employ a linear basis for the annihilator, which will 
form a linear systems of partial differential equations of finite type. The outline of this 
method is reasonably clear, but full details remain to be worked out. In particular, the 
existence of corresponding bounds on the dimensions of invariant subspaces for nonlinear 
partial differential operators is not known. 

6. Applications to Differential Equations. 

A principal application of our theory is to find explicit solutions to both linear and 
nonlinear differential equations. In the theory of quasi-exact solvability, one considers an 
eigenvalue problem 



which, in physical applications, is the stationary Schrodinger equation. A linear differential 
operator Q is said to be quasi-exactly solvable if it leaves a finite-dimensional subspace of 



zK^[u] = D^F[z,u] 



(5.12) 



Q[u] = Alt, 



(6.1) 
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wave functions M. <Z invariant, so that Q[u] e V{M.). In this case, restricting to the 
subspace Ai = • • • , /^.(a;)}, whereby 

u{x) = c^fiix) H c^frix), (6.2) 

reduces the eigenvalue problem (6.1) to a linear eigenvalue problem 

Ac= Xc, for c = (c;^, . . . , c^)"^. (6-3) 

In this way, one reduces the solution of the differential equation (6.1) to an algebraic 
problem. More generally, one can consider nonlinear differential operators Q[u\ e Q{M.) 
that leave M invariant, in which case the restriction of (6.1) to the subspace leads to a 
system of nonlinear algebraic equations 

H(c) = Ac (6.4) 

for the coefficients in (6.2). 

For evolution equations, the basic idea underlying the method of nonlinear separation 
goes back to Galaktionov and his collaborators, [11-15], King, [26], and Fushchych and 
Zhdanov, [9, 10]. The basic idea appears in the following theorem. 

Theorem 6.1. Consider an evolution equation 

Ut = Q[u]. (6.5) 

Suppose the right hand side Q[u\ e Q{M.) preserves a Gnite-dimensional subspace Ai, that 
is 

r I 

Q[u] = HiiL.H • • • , LM) fi + Yl ■ KH 
i=l 1^=1 

where the e C'^(M^) are arbitrary analytic functions, and where the e ^(-^) are 
arbitrary operators. Then there exist "separable solutions" of the evolution equation taking 
the form 

n 

u{x,t) = J2'Piit)fiix), (6.6) 

i=l 

if an only if the coefficients (p^, . . . ,(fi^ are solutions to the dynamical system 

dep. 



dt 



Hi{Vi.---.Vn)^ z=l, ...,n. (6.7) 



The proof is immediate from Theorem 5.3. More generally, one can replace the evo- 
lution equation (6.5) by a time-dependent Schrodinger equation 

iu^ = Q[u], 

resulting a complex first order dynamical system, or a wave-type equation 

Utt = Qbj], 
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which reduces to a second order dynamical system, with dPip^/dt^ on the left hand side of 
(6.7). Indeed, one can apply the method to dynamical partial differential equations of the 
general form 

T[u]^Q[u], (6.8) 

where T is a linear ordinary differential operator in t, or, even more generally, a time- 
dependent linear combination of operators Tf,[d^u/ dt^) where each e Q{M.) leaves M. 
invariant. An interesting class of examples are the equations of Fuchsian type, studied in 
depth by Kichenassamy, [24] in connection with blow-up and Painleve expansions, [25], 
in which T[u] is a constant coefficient polynomial in the scaling operator S = tdf. 

Remark: See Cherniha, [5], for an even more general nonlinear separation ansatz in 
which the basis functions fi{x, t) (and hence the module) are also allowed to depend on t. 

Example 6.2. We can already give many examples of non-linear QES evolution 
equations in two space variables by using the generators (4.8). Consider for example the 
simplicial subspace 

The dual basis of operators is given in (4.9). Therefore, the most general second order 7^ 
invariant evolution equation takes the form 

u^ = A + Bx + Cy + Dx'^ + Exy + Fy^, 
where A, S, C, D, F are arbitrary (linear or nonlinear) functions of the dual operators 

^xxi '^xy ^yyi "^^xx V^xyi 

Uy - xu^y - yuyy, \x^u^^ + xyu^y + \y'^Uyy - xu^ - yuy + u. 

Higher order evolution equations are obtained by adding in arbitrary annihilators. Re- 
placing by u^^ or other types of linear temporal differential operators, leads to wave and 
more general types of equations that leave the indicated subspace invariant. It follows, for 
example that the evolution equation 

«t = i a^^^L - ^y^xy - y\y + 1 y^^ly 

admits the solutions given by 

uix, y, t) = + k^e'* xy + y^ + h + + hy^ 

where the are arbitrary constants determined by the initial conditions. This is a very 
simple example which was chosen such that the dynamical system governing the Cj(t) was 
decoupled and therefore integrable by quadratures. 

Example 6.3. As a more substantial example, consider the following rotationally 
invariant evolution equation 

u^ = {A- B)[V{u'^Vu) - Au'^+^] + Cu + d, u = u{x,t), (7=1,2 (6.9) 
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in n space variables x = {x^, . . a;"). We claim that the subspace 

is invariant for a suitable choice of a^^, To see this let us define 

R[u\ = grad(«'^ gradw) — Au'^'^^, 
and note the following identities: 

R[u + v]= R[u] + R[v] + \ / ^ M o 

1^ [A — A)[uv{u + v)\, (7 = 2. 

Also let us note that for a e M", 

R[e^-^] = [(a + l^aW^ -A] e^'^ 
Hence, taking a-^, 0.2 so that llajp = ^/((7+ 1), z = 1,2, we will have 

i?[C]^tt]^] = R[c2U2\ = 0, for all C]^,C2 € 

Of course i? does not annihilate all oi M.\ we are left with quadratic cross-terms. Indeed, 
taking u = c^u^ + C2M2 we have for cr = 1 

R[u\ = R[CiUi] + R[C2U2] + CiC2{A — A)[UiU2]. 

Hence, taking B — \\di + CK2IP ^^^^ ensure that 

(A - B)[R[u]] = 0. 

For cr = 2, we need to take B = \\2a-^ + a2|p, which is the same as + 2a2|p. Indeed, 
we have: 

R[u] = c-^C2iA-A)[u^U2iu^+U2)] = C^{e^^^'+^^^-^ + e^^'+^^^>''), 

and the right hand side is manifestly annihilated hy A — B. It isn't hard to check that 
R[u + C3] differs from R[u] by a linear combination of Ui, U2, and a constant. Hence, 
{A — B)[R[u + c^]] continues to lie in M., and therefore the right hand side of equation (6.9) 
preserves A4. 

We could also consider the extended module 

M = {uq = 1, = e°l■^ U2 = e°2-^, w_i = e""^-^, u_2 = e-^^'^}. 

The same reasoning as above can be applied to show that R[u], where w is a general 
element of Ai, is the sum of a linear combination of i G {1, 2, —1, —2} and a number 

of cross-terms of the form constant times u-Uj, i ^ j, i,j e {0, 1, 2, —1, —2}. We therefore 
require that A — B annihilate u^U2i u^u_2, and their reciprocals. For cr = 1, this constraint 
requires 

S = 2||a]^||^ = 2||a]^||^, and • ^2 = 

so that and a2 are perpendicular. For cr = 2, the conditions are: 

5 = 5||a]^|p = 5||a]^|p, and ■ a2 = 

and the conclusions are identical. We remark that, by making use of the rotational invari- 
ance of Q, we can without loss of generality assume that and that a2 - x = x . 
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7. Conclusions. 



In this paper, we have developed a general theory of invariant subspaces, and shown 
how these results can apply to provide nonlinear separation of variables ansatze for a 
wide variety of linear and nonlinear partial differential equations. These operators further 
provide a nonlinear generalization of the quasi-exactly solvable operators of importance in 
the algebraic approach to quantum mechanical systems. A number of interesting problems 
warrant further development of this method. 

(1) A significant mystery is the connection of this method with the Lie algebraic approach 

to quasi-cxactly solvable modules. For example, the linear operators preserving 
the simplicial subspacc (4.6) lie in the universal enveloping algebra of a standard 
realization of the Lie algebra gl(n, M), cf. [16,46]. However, this is not evident 
from the form of the affine annihilator and annihilator. Thus, a fundamental 
issue is which subspaces admit such a Lie algebraic interpretation of the space of 
differential operators that leave them invariant. 

(2) The "inverse problem" of characterizing the invariant finitc-dimcnsional subspaces 

for a given linear or nonlinear operator is of critical importance. The extension of 
Svirshchevskii's symmetry approach, [49], is under investigation. 

(3) Applications of our analytical Wronskian methods to the algebraic context of D- 

module theory, [2, 31, 32], looks quite promising. In particular, the characteriza- 
tion of the analytic or polynomial annihilators of subspaces of rational functions 
would be a particularly interesting case. 

(4) The formulae for the affine annihilators and annihilators are often extremely compli- 

cated, even for relatively simple subspaces. (As an example, the reader is invited 
to write down the formulae for = {x^, xy, y^}-) Moreover, in the final formulae 
(5.8), (5.7) for the linear and nonlinear operators leaving the subspace invariant, 
one may encounter a significant amount of simplification and lowering of order. 
Thus, it would be important to characterize low order operators in V{A4) and 
Q{M), as well as "simple" operators of physically relevant type, e.g., elliptic, 
constant coefficient, Lorentz invariant, etc. 

(5) Extensions of these methods to finite difference operators, building on the work of 

Turbiner, [46], can be profitably pursued. 
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